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^vq Abstract 

i— —i In this article we study convex integer maximization problems with com- 

posite objective functions of the form f(Wx), where / is a convex function 
on W 1 and W is a d x n matrix with small or binary entries, over finite sets 

^ S C Z n of integer points presented by an oracle or by linear inequalities. 

Continuing the line of research advanced by Uri Rothblum and his col- 
leagues on edge-directions, we introduce here the notion of edge complexity 
of S, and use it to establish polynomial and constant upper bounds on the 

,_l number of vertices of the projection coiw(WS) and on the number of linear 

J> optimization counterparts needed to solve the above convex problem. 

Two typical consequences are the following. First, for any d, there is a 

\q constant m{d) such that the maximum number of vertices of the projection of 

any matroid S C {0, l} n by any binary d x n matrix W is m{d) regardless of 
n and S; and the convex matroid problem reduces to m(d) greedily solvable 
linear counterparts. In particular, m(2) = 8. Second, for any d, I, m, there is a 
constant t(d; I, m) such that the maximum number of vertices of the projection 
of any three-index I x mx n transportation polytope for any n by any binary 
dx(lxmxn) matrix W is t(d; I, m); and the convex three- index transportation 

rN problem reduces to t(d; I, m) linear counterparts solvable in polynomial time. 
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1 Introduction 

In this article we study convex integer maximization problems and the closely related 
projections of the sets of feasible points. Let S C 7L n be a finite set of integer 
points, let conv(S') C M. n be its convex hull, let W be a d x n integer matrix, and 
let / : M d — > R be a convex function. We study the problem of maximizing the 
composite function f(Wx) over S and the projection of conv(S') by W into Mr, 
namely, 

max {f(Wx) : x G S} (1) 

and 

conv(W / 5') = coiav{Wx : x £ S} = {Wx : x G conv(S)} C R d . (2) 

The sets S we consider arise in two natural contexts. First, in combinatorial 
optimization, in which case S C {0, l} n has some combinatorial structure and might 
be presented by a suitable oracle. Second, in integer programming, where 

S := {x G Z n , Ax = b , I < x < u} (3) 

is the set of integer points satisfying a given (standard) system of linear inequalities. 

The optimization problem can also be interpreted as a problem of multicri- 
teria optimization, where each row of W gives a linear criterion W{X and / compro- 
mises these criteria. We therefore call W the criteria matrix or weight matrix. 

The projection polytope conv(WS) in ^ and its vertices play a central role in 
solving problem ([T]): for any convex function / there is an optimal solution x whose 
projection y := Wx is a vertex of conv(WS). In particular, the enumeration of all 
vertices of conv(WS) enables to compute the optimal objective value for any given 
convex function / by picking that vertex attaining the best value f(y) = f{Wx). 
So it suffices, and will be assumed throughout, that / is presented by a comparison 
oracle that, queried on vectors y, z G W 1 , asserts whether or not f(y) < f(z). 

This line of research was advanced by Uri Rothblum, to whom we dedicate this 
article, and his colleagues, in several papers including [21 El H5] . and culminated in 
the edge-direction framework of |16j . see also [151 Chapter 2]. In this article we 
continue this line of investigation, and take a closer look on coarse criteria matrices; 
that is, we assume that the entries of W are small, presented in unary, or even 
bounded by a constant and lie in {0, 1, . . . ,p}. In multicriteria combinatorial opti- 
mization, this corresponds to the weight Wy attributed to element j of the ground 
set {1, . . . , n} under criterion % being a small or even {0, 1} value for all i,j. 

Here is a typical result we obtain in convex combinatorial optimization, where 
S C {0, 1}™ is the set of indicating vectors of bases of a matroid over {1, . . . , n}. 
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Corollary 4.7| The maximum number of vertices of conv(WS) for any n, any 

matroid S C {0, l} n and any binary 2xn matrix W, is 8; and max{f(Wx) : x G S} 
can be solved for any convex f by greedily solving 8 linear counterparts over S . 

Note that, in contrast, there are sets S C {0, l} n and {0, l}-valued 2xn matrices 



W such that the projection conv(WS) has D,(-^n) vertices, see Theorem 3.1 

More generally, we show in Section 4 that for any d and p there is a constant 
m(d,p) such that the maximum number of vertices of the projection of any matroid 
S C {0, 1}" by any {0, 1, . . . ,p}- valued dx n criteria matrix W is m(d,p) regardless 
of n and S. Thus, m(2, 1) = 8 but the precise values of m(d,p) for all higher d and 
p are unknown and their determination poses a challenging open problem. 

Moving to convex integer programming, we have the following typical result 
concerning multi-index transportation problems. 



Corollary 5.8 For every fixed d } p,l } m, there exists a constant t(d,p;l,m) such 
that, for every n, every integer line-sums a^k, cy, and every {0, 1, . . . ,p}-valued 
d x (7 x m x n) matrix W , the projection cortv{WS) of the set of I x m x n tables 

S . \ X G TLj^ . ^ •Kj,j,k ®j,k ) ^ . 3*1, j,k j ^ . ^i,j,h \ j (4) 

L i j k ) 

has at most t(d,p; I, m) vertices regardless ofn and the line sums. Moreover, all these 
vertices can be enumerated, and for any convex f the problem m&x{f(Wx) : x G S} 
solved, by solving t(d,p; I, m) linear counterparts over S, and in polynomial time. 

Other applications, to convex totally unimodular integer programs, network flows 
and more generally transshipment problems with convex multicriteria objective func- 
tions, and vector partition problems, are discussed in Section 5. 

These results and applications follow from a general result which we now describe. 
Let S C 7h n be a finite set of integer points. We define the edge complexity of S to be 
the smallest nonnegative integer e(S) such that every edge of the polytope conv(S') 
is parallel to some vector v G Z n with ||i>||i < e(S). For a matrix W we use the 
notation H^Hoo := maxy |Wij|. We establish the following theorem. 

Theorem 2.3 Fix any d. Let S C Z" be a finite set of integer points and W an in- 
teger dxn matrix. Then conv(WS) has polynomially many O ((e(S) ■ || W^Hoo)^ 
vertices. Moreover, if S is presented by a linear- optimization oracle, and endowed 
with an upper bound e on e(S), then all vertices of comr(WS) can be enumerated, 
and meix{f(Wx) : x G S} solved for every convex f : M. d — > K, in polynomial time. 



The dependency in Theorem 2.3 on the unary size of e(S) and W cannot be 
relaxed: in Theorems 3^2 and 3^ we show, respectively, that already projections 
into the plane M 2 can have exponentially many vertices, even when S C {0, l} n and 
hence the edge complexity of S satisfies e(S) < n, or when W is a {0, l}-matrix. 
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Theorem 2.3 continues the line of research taken in [21 [31 El [91 [131 HB] implicitly or 
explicitly. But in these papers, a complete set of edge-directions of conv(S') in M. n was 
required as part of the input, whereas here, only a bound on the edge complexity is 
needed. Consequently we get constant bounds on the number of vertices and number 



of linear counterparts needed, in Corollaries 4/7 and ^8 and other applications, even 
in situations where conv(S') has exponentially many edge-directions. 

The rest of the paper is organized as follows. In Section 2 we discuss edge com- 



plexity and projections of polytopes, prove our general Theorem 2J3 which leads to 
our polynomial and constant bounds on the number of vertices of projections, and 
describe the algorithm underlying it. In section 3 we construct polytopes whose 
projections have large number of vertices, in contrast with Theorem 2.3 In partic- 



ular, we show in Theorems 3^ and 3.3| respectively, that projections into the plane 
of {0, l}-sets S or by binary matrices W can have exponentially many vertices. 

In Section 4 we focus on matroids in convex combinatorial optimization. We 



show in Theorem 4.1 the existence of bounds m(d,p) on the number of vertices of 
the projection by any {0, 1, . . . ,p} valued d x n matrix W of any matroid S over a 
ground set of any size n, and reduce the convex matroid optimization problem to 
m(d,p) linear counterparts. We then restrict attention to binary matrices W and 
give sharper upper and lower bounds on m(d) := m(d, 1) in Theorem 4.6 For d = 2 
these coincide and we conclude Corollary ~AJ\ with m(2) = 8. We conclude Section 



4 with a detailed discussion of bicriteria matroid optimization in Example 4.8 

In Section 5 we proceed to convex integer programming. We discuss Graver 
bases and their connection to edge complexity. In Section 5.1 we prove Theo- 
rem |5.2| on convex totally unimodular integer programs, and discuss applications 



to convex multicriteria transshipment in Corollary |5.3| and vector partition prob- 
lems in Corollary |5.4 In Section 5.2 we prove Theorem 5.7 on convex n-fold integer 
programming, discuss the application to muti-index transportation problems and 
demonstrate Corollary 5.8, and prove an upper bound on the maximum number 



t(d; /, m) := t(d, 1; /, m) of vertices in any projection of 3-way tables by binary crite- 



ria matrices in Proposition 5.9 We conclude in Section 6 with some open problems. 



Before proceeding we set some computational complexity notation. The binary 
lengthof an integer z is the number (z) = 0(log(|z|+l)) of bits in its binary encoding. 
The binary length (x) of an integer vector x is the sum of binary lengths of its entries, 
and the length (E) of a set E of integer vectors is the sum of lengths of its elements. 
Our algorithms have running times which are typically polynomial in n and ||W||oo; 
in the binary length of the rest of the data if any, which in integer programming 
is typically the binary length (b, I, u) of the right-hand side and lower and upper 
bounds; and they perform a polynomial number of queries to the comparison oracle 
presenting / and the linear-optimization oracle presenting S when relevant. 



5 

2 Edge complexity and projections of polytopes 



We begin with some preliminaries on edge-directions and zonotopes. A direction 
of an edge (1-dimensional face) e of a polytope P is any nonzero scalar multiple of 
u — v with u, v the vertices of e. A set of all edge- directions of P is a set which 
contains some direction of each edge of P. The normal cone of a polytope P C R n 
at its face F is the relatively open cone of those linear functions maximized 
over P precisely at points of F. A polytope Z is a refinement of a polytope P if 
the normal cone of every vertex of Z is contained in the normal cone of some vertex 
of P. The zonotope generated by a set of vectors E = {ei, . . . ,e m } in R d is the 
following polytope, which is the projection by E of the cube [—1, l] m into R d , 

'^2\ i e i : Xi = ±1 > = ^2 [- e i> e i] C R d . 

i=l J i=l 

We have the following two lemmas, see for instance [8] or [151 Chapter 2] and the 
references therein for more details. The lemmas are illustrated in Figure [T] below. 

Lemma 2.1 Let P be a polytope and let E be a finite set of all edge-directions of 
P. Then the zonotope Z := zone(E) generated by E is a refinement of P. 

Lemma 2.2 The number of vertices of any zonotope Z := zone(E) generated by a 
set E of m vectors in R d is at most 2 Ylt=o ("V*) ■ For every fixed d, there is an 
algorithm that, given E C Z d , outputs every vertex u of Z := zone(E) along with 
some h u G 1/ maximized over Z uniquely at u, in time polynomial in m and (E). 

We next introduce a new concept which is central to our development. Let the 
edge complexity of a rational polytope P be the smallest nonnegative integer e(P) 
such that every edge of P is parallel to an integer vector v with ||i>||i < e(P). Let 
the edge complexity of a finite set S of rational points be e(S) := e(conv(S')). For 
matrix W let H^Hoo := maxy A linear- optimization oracle for S is one that, 

queried onw£ Z n , solves the linear optimization problem m&x{wx : x G S}. 

We are now ready to prove our theorem on edge complexity. Its algorithmic 
statement is furnished by the algorithm outlined following the proof of the theorem. 

Theorem 2.3 Fix any d. Let S dTI 1 be a finite set of integer points and W an in- 
teger dxn matrix. Then conv(WS) has polynomially many O ((e(S) ■ || W / || 00 ) a! ^ _1 ^) 
vertices. Moreover, if S is presented by a linear- optimization oracle, and endowed 
with an upper bound e on e(S), then all vertices of comr(WS) can be enumerated, 
and meix{f(Wx) : x G S} solved for every convex f : R d — > R, in polynomial time. 
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Figure 1: A zonotope refining a polytope and unique maximizers at their vertices 

Proof. Let e := e(S) and q := e • ||I^||oo- Let P := comr(WS) be the projection 
polytope. Now, consider any edge / = [u, v] of P. Then there must be an edge g of 
conv(S') whose projection Wg is not a singleton and satisfies Wg C /. Since g is an 
edge of conv(S'), it is parallel to some feeZ" with \\h\\i < e(S) < e. Therefore Wh 
is an edge-direction of / and satisfies Wh G {0, ±1, ... , ±q} d , since, for % — 1, . . . , d, 

n n 

\Wih\ < ^{WijWhj] < H^IU^I^-I = H^llooll^lli < q. 



Since this is true for every edge of P, it follows that E := {0, ±1, ... , ±q} d is a 
set of all edge-directions of P. Let Z := zoiae(E) be the zonotope generated by 



the (2q + 1) -set E. Since Z is a refinement of P by Lemma 2.1, the closure of 



each normal cone of P is the union of closures of normal cones of Z. Therefore the 



number of vertices of P is bounded by that of Z, which by Lemma is at most 

d-l 



k=0 ^ ' k=0 ^ 



d(d-l)^ 
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We continue with the algorithmic part of the theorem (see also the description of 
the algorithm following this proof). Let e be the given upper bound, let again 



q := e ■ HWUoo and E := {0, ±1, . . . , ±q} d . By Lemma 2.2 in polynomial time we 
can produce every vertex u of the zonotope Z := zone(-E) along with some h u G Z d 
maximized over Z uniquely at u. Now consider any vertex v of P := comr{W S). 
Since Z refines P by Lemma [2. 1[ there is a vertex u of Z such that the normal cone 
of Z at u is contained in the normal cone of P at v. Therefore, h u is maximized 
over P uniquely at v. Let x u G S 1 be a maximizer of /i^W 7 G Z n over S 1 . Then, 
since for every x G conv(S') and its image y := Wx G convfWS 1 ) = P we have 
^u2/ = h u Wx, it must be that u = VTa; u . So the algorithm proceeds as follows (see 
description of the algorithm following this proof). For each vertex u of Z, use the 
linear-optimization oracle of S to obtain a maximizer x u G S of /i^W G over S, 
and collect all such x u in a set X C Z n and all such v = Wx u in a set 1/ C Z^. Then 
V is the set of vertices of P. Moreover, an optimal solution x* is any point in X 
attaining maximum f(Wx), which can be found using the comparison oracle of /. □ 



The following furnishes the algorithm used and analyzed in the proof above. 



Projection and Convex Maximization Algorithm (of Theorem 2.3) 



Input: Finite set SCZ™ presented by a linear-optimization oracle, 
upper bound e on its edge complexity e(5*), integer d x n matrix W, 
and convex function / : M. d — > K presented by a comparison oracle. 



1. Preprocessing: 

Let q := e • H^Hoo- Construct the zonotope Z = zone{0, ±1, . . . , ±q} d , 
and for each vertex u of Z, some h u G Z d maximized over Z uniquely at u. 

2. Projection: 

Let X := 0, V := 0. For each vertex u of Z, obtain x u G S maximizing 
the linear function h u W over 5* and set X := X U {x u }, V := V U {Wx u }. 

3. Optimization: 

Obtain x* G X attaining maximum value of f{Wx) over all i6l. 

Output: The set V of vertices of the projection conv(WS) 
and the optimal solution x* of ma.x{f(Wx) : x G S*}. 
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3 Projections with exponentially many vertices 

In this section we construct sets S of integer points and weight matrices W such 
that the projections conv (WS), already into the plane IR 2 , have many vertices, 
contrasting our polynomial and constant bounds of the previous and next sections. 

We begin by considering the case with both S C {0,l} n and W G {0, l{ 2xn 
binary valued. Then the projection satisfies WS C {0,1,..., n} 2 and hence the 
number of vertices is 0(n 2 ) and cannot be exponential. We show it can be Q (y/ri). 

Theorem 3.1 For every n there exist set S C {0, l} n and matrix W G {0, l{ 2xn 
such that the projection ofconv(S) by W into the plane has more than \ \n vertices. 



Proof. Let k := y |n so that n > k + k 2 . Define a matrix W with each of its first 

k columns being the unit vector l x G IR 2 and each of its next k 2 columns being the 
unit vector I2 G IR 2 . Define a set S in terms of the unit vectors lj G IR n by 

s ■■= |e^+ e h ■ o<i<k 

^ j<i k+l<j<k+i 2 

Then WS = {(i,i 2 ) : < i < k} is a set of k + 1 > \J\ n points in convex posi- 
tion, which comprise the vertices of the projection of conv(S') by W into the plane. □ 

Next we construct several sets of integer points in polytopes whose projections 
into the plane IR 2 have exponentially many vertices, in contrast with Theorem 2.3 



First, we construct integer programs defined by {—1,0, 1} matrices, and weight 
matrices W G Z 2xn with binary size (W) polynomial in n, such that the integer 
hulls are {0, l}-polytopes, and their projections have exponentially many vertices. 
Note that the edge complexity of any {0, l}-polytope is at most n, since any edge is 
parallel to the difference between its two vertices which is a {—1, 0, 1 {-valued vector. 



Theorem 3.2 For every k there are m,n < 2k 2 and matrices A G {0,±l{ mxri and 
W G l}+ n with log llWlloo < 2k, such that the projection by W of the {0, \}-polytope 
Q := conv{a; G Z, n : Ax = 0,0<x<l{ into the plane M. 2 has exactly 2 k vertices. 

Proof. Consider the following system of m := 3( 2 ) equations in n := k + 4(„) binary 
variables Xj, < i < k and Xij, Ui t j,Vi t j, Zij, < i < j < k, 
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and let A be the corresponding mxn matrix. Let W be the matrix whose rows W\, W2 
give the following linear functions in the vector of variables x := (a^, x it j, u^j, i>y, Zy), 

k-1 k-1 

Wix := ^2'xi, w 2 x := ^4^ + 2 V} 2 i+j Xij . 

i=0 i=0 0<i<j<k 

Clearly m,n < 2k 2 , A is {0, ±l}-valued, and (|| W||oo) < 2k, as claimed. 

It is easy to see that the equations force that, for any < i < j < k, if Xi = or 
Xj = then = 0, whereas if Xj = Xj = 1 then xy = 1. Thus, every choice of the 
partial variable vector x := (xo,xi, . . . ,Xk-i) G {0, l} fc can be uniquely extended to 
a vector x G {0, l} n satisfying the above system of equations. So the integer points 
in Q are in bijection with such choices x and hence Q has 2 k vertices. Moreover, for 
each such choice x, the equations force x^j = XiXj for all i,j, and hence 



k-1 \ " fc-1 

(wix) 2 = [^2 i x i I = ^(2 i ) 2 x 2 + 2 ^ 2 4 x 4 2^ 

\i=0 / i=0 0<i<j<k 

k-1 

i=0 0<i<j<k 

When x runs over {0, l} fc , the value W\x runs over all values 0, 1, . . . , 2 k — 1. So the 
set of integer points of Q projects to the set {(£, t 2 ) : t = 0, 1, . . . , 2 k — 1} of 2 k 
points in convex position, comprising the vertices of their convex hull conv(W / S). □ 

The last theorem of this section constructs integer programs defined by {—1, 0, 1, 2}- 
valued constraint matrices, and {0,l}-valued 2 x n weight matrices W, such that 
the number of vertices of the projections of the integer hulls is exponential in n. 

Theorem 3.3 For every k there are m,n < 4k 4 , matrices A e { — 1, 0, 1, 2} mXTl and 
W G {0, \} 2xn , and lower and upper bound vectors I G {— oo, 0} n , u G {1, oo} n , such 
that the integer hull Q := conv{x G Z n : Ax = 0,l<x<u}isa polytope whose 
projection by the binary matrix W into the plane M. 2 has exactly 2 k vertices. 



Proof. As in the proof of Theorem 3.2[ consider again the following system of 



equations in binary variables X®, < i < k and x®j, Uij, Vy, Zij, < i < j < k, 

x °i ~ x lj ~ u i,j = > x °j ~ x lj ~ v i,j = ^ x °i + x °j ~ x lj ~ z i,j = > 0<i<j<k. 

For each variable x° introduce 2i additional variables x\ without bounds and 2i 
additional equations 2x\~ x — x\ — 0, r = 1, . . . , 2i. Likewise, for each variable X® • 
introduce i + j + 1 additional variables x\^ without bounds and i + j + 1 additional 
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equations 2x T ~^ — x\j = 0, r = 1, . . . , i + j + 1. Then there are n = variables and 
m = equations with {0, ±1, 2}-valued defining matrix A. 

Let W be the {0, l}-valued matrix whose rows w±, W2 give the linear functions 



fc-i fc-i 



WlX := J^x-, w 2 x := ^ rc * +J+ ' 



i=Q i=0 0<i<j<k 

It is easy to see that the new equations force on the new variables the relations 



x[ = 2 r x i and x[ ■ = 2 r x\ • for all i,j,r. Therefore, we obtain 



fc-i fc-i 

» ,3 



Wl x = ^2^°, w 2 x := J^4X° + 2 ^ 2 i+J x° 



0<i<j<fc 



The proof now proceeds as that of Theorem that is, every choice of the partial 
variable vector x := (x®, . . . , G {0, l} fe can be uniquely extended to a vector 
x E Z" satisfying the system of equations, that is, an integer point in Q, and these 
points project down onto the set {(t, t 2 ) : t = 0, 1, . . . , 2 k — 1} of 2 k points in convex 
position, which comprise the vertices of their convex hull conv(WS). □ 



4 Convex matroid optimization 

In this section we study the convex multicriteria optimization problem ([!]) and the 
corresponding projection ^ over matroids. We identify a matroid with its set 
of bases, so we call S C {0, 1}™ a matroid if it is the set of (indicating vectors 
of) bases of a matroid over {1, . . . , n}. (We restrict attention to bases, but similar 
results hold for independent sets.) In particular, our results below apply to spanning 
forests in graphs. An independence oracle for S is one that, queried on vector 
x G {0, l} n , asserts whether or not x is independent in S, that is, whether or not 
supp(x) C supp(2;) for some z G S. As is well known, when a matroid S is presented 
by an independence oracle, the linear optimization problem max{u>x : x G S} can 
be easily solved by the greedy algorithm. See [IB] for more details on matroids. 



As shown in Theorem 3.1, when projecting a binary set S C {0, l} n by a binary 
matrix W G {0, l} dxn ; even to the plane d = 2, the number of vertices can grow with 
n. In contrast, we now show that for matroids, the number of vertices is constant. 

Theorem 4.1 For every d and p, the maximum number of vertices of conv (WS) 
for any n, any matroid S C {0, l} n , and any {0, 1, . . . ,p}-valued dx n matrix W , is 
a constant m(d,p) which is independent of n, S and W. Moreover, if S is presented 
by an independence oracle then m&x{f(WX) : x G S} can be solved for any convex 
f by greedily solving m(d,p) linear counterparts over S, and in polynomial time. 
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Proof. Let S be any matroid. Then for every edge e = [u, v] of the matroid polytope 
conv(S') we have that u — v = lj — lj is equal to the difference of some two unit 
vectors in IR n ; here is the proof, included for completeness. Let B u := supp(u) and 
B v := supp(f) be the corresponding bases. If B u \ B v = {i} is a singleton then 
B v \ B u = {j} is a singleton as well in which case we are done. Suppose then, 
indirectly, that this is not the case. Let h G IR n be uniquely maximized over P at e, 
and pick an element i in the symmetric difference B U AB V := (B u \ B v ) U (B v \ B u ) 
of minimum value hi. Without loss of generality assume i G B u \ B v . Then there is 
a j G B v \ B u such that B := B u \ {i} U {j} is also a basis. Let x := 1b G S be the 
indicator of B. Now \B U AB V \ > 2 implies that x ^ u and x ^ v. By the choice of 2 
we have hx = hu — hi + /ij > hu. So x is also a maximizer of h over P and hence 
x G e. But no {0, l}-vector is a convex combination of others, a contradiction. 

Thus, every edge of conv(S') is parallel to some lj — lj and hence S has constant 
edge complexity e(5*) = 2. Also, the criteria matrix satisfies H^Hoc — V- Finally, a 
linear-optimization oracle for S is realized from the independence oracle presenting 



it by the greedy algorithm. The theorem now follows from Theorem 2.3 □ 



As mentioned in the introduction, the determination of the constants m(d,p) is 
a challenging open problem. We now restrict attention to binary criteria matrices, 
that is, {0, l}-valued d x n matrices W. Let m(d) := m(d, 1) be the corresponding 
maximum number of vertices of the projection conv(WS) of any matroid S by any 
binary matrix. Introduce the zonotope Z(d) := zone{ — 1, 0, l} d and let z(d) be the 
number of vertices of Z(d). We next provide an upper bound on z(d) and m(d). 

Lemma 4.2 The maximum number of vertices of the projection conv(iyS') into W 1 
of any matroid S by any binary d x n matrix W satisfies the upper bound 

d-l 



m{d) < z{d) < 2j2( U3d k 3) ) 

k=0 ^ ' 



Proof. The proof of Theorem 2.3 implies that every edge of conv(WS) is parallel 



to W h for some edge-direction h of conv(S'). Since S is a matroid, the proof of 



Theorem 4.1 implies that every edge of conv(5') is parallel to some 1, — lj. Since 
W G {0, l} dx?1 , it follows that a set of all edge-directions of conv(W / S') is given by 

{Wh : h = li-lj, 1 < i < j < n} C {-1,0, l} d . 

So Z(d) = zone{ — 1, 0, l} d refines coiav{WS) and the latter has at most m(d) < z(d) 
vertices. Now, Z(d) is homothetic to the zonotope generated by the following set 
containing only one of each antipodal pair of nonzero vectors in { — 1, 0, l} d , 



{v G { — 1, 0, l} d : v 7^ and the first nonzero entry of v is 1} 
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Since this set has \{?> d — 1) elements, the bound now follows from Lemma 



2.2 



We proceed to develop some lower bounds on m(d) (and hence also on z(d)). 
The uniform matroid of rank r and order n is the matroid having all r-subsets of 
{1, . . . , n) as bases. We identify it with the set 5" C {0, l} n of all vectors x G {0, l} n 
with |supp(x)| = r. Let u(d) be the maximum number of vertices of the projection 
convfWS"™) of any uniform matroid under any binary d x n matrix W. Clearly, 
u(d) < m(d) for all d. For positive integers k, d, let Wjj be the d x k2 d binary criteria 
matrix whose k2 d columns consist of k copies of each vector in {0, l} d . Finally, let 
u^(d) be the number of vertices of coiav{W^S^ 2 ); note that it is well defined and 
independent of the order of columns of W% by the symmetry of uniform matroids. 
Clearly, u^(d) < u(d) < m(d) for all k, r, d. Beyond providing lower bounds on m(d), 
the determination of u(d) and u^(d) is an interesting problem in its own right; the 
following three statements characterize the vertices of conv(W^S^ 2d ) and determine 
Ur(d), respectively, for k > r, k = 1 and r = 2, and k > 2 and r = k + 1, for all d. 

Proposition 4.3 Fork > r, conv(W^S^ 2d ) = [0, r} d is a cube and hence u^(d) = 2 d . 

Proof. Clearly the image of S^ 2 under W% satisfies W^S* 2 C {0, 1, . . . ,r} d and 
hence is contained in the cube. On the other hand, every vertex v of the cube [0, r] d 
is a multiple v = r ■ u for some u G {0, l} ra ; letting J be any set of r < k distinct 
indices of columns of W§ which are equal to u, and letting x := Z~2jej lj e $r 2 > we 
find that v = W^x G W^S l f' d is in the image. So the image is contained in [0,r] d 
and contains all its vertices and hence conv(W^S^ 2d ) = [0,r] d as claimed. □ 

We assume next k < r and begin with k — 1 and r = 2. Let V := {0, 1, 2} d and 
for % — 0, 1, . . . , d let Vi be the set of vectors in V with precisely i entries equal to 1. 

Proposition 4.4 For k = 1 and r = 2, the image satisfies W^S^ = V\ Vq, the set 

of vertices of its convex hull conv(W / }5'| ) is equal to V\, and therefore u\(d) = d2 d ~ 1 . 

Proof. Let n := 2 . The image W^S^ consists of all sums u + v of distinct pairs 
u,d£ {0, l} d . Each such sum u+v is in V and for some i we have Ui ^ Vi so Ui+v j = 1 
and hence it + 1> G" Vo- Conversely, consider z G V \ V with ^ = 1« + 2^jeJ ^-^i- 
Then z = u + v with it := J^ie/uJ •"•* an< ^ 11 := J^jeJ distinct since 7^0. 

We proceed with the claim about vertices. Consider any v G W^S^Vi. By what 
we just proved, v G" V and so Vi = 1 for some i. Since i; G" V\ we have that d + lj 
and u — lj are not in Vb hence in the image W\SV[ , and v = \ ((v + lj) + (v — lj)) 
is a convex combination of them and hence not a vertex. Conversely, consider any 
v G V\. Then i>j = 1 in precisely one entry. Define a vector h G Z d by 

— 1, if = 0, 
/lj := ^ 0, if Vj = 1, 
1, if vj = 2. 
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Then hu < hv for all u G V, with equality if and only if u G {v — lj , v , v + lj}. 
But v — lj , v + lj G Vo so v is the unique maximizer of /i over W^JS^ hence a vertex. 
So the number of vertices of coiaviW^S^) is w^cf) = |Vi| = d2 d_1 as claimed. □ 

We proceed with the case of r = k + 1 for any k > 2. Let V := {0, 1, . . . , r} d and 
for i = 0, 1, . . . , d let Vi be the set of vectors in V with i entries not equal to of r. 

Lemma 4.5 For k > 2 and r = A; + 1, the image satisfies W^S^ 2d = V\Vq, the set 
of vertices of comr{W^S^ 2 ) consists of those vectors v G V\ which have one entry 
equal to 1 or r — 1 and all other entries equal to or r, and therefore u^(d) = d2 d . 

Proof. Let n := k2 d . The image W^S™ consists of all sums z = Y7i=i ut °f r vectors 
u l G {0, l} d with at least two of them distinct. Clearly each such sum z is in V and 
for some i, j, h we have u l h ^ u 3 h so u\ + u J h = 1 and hence 1 < < r — 1 so z £ Vq. 
Conversely, consider z G V \ V . Then z = Yll=i 2~2j£ii f° r some pairwise disjoint 
Ii with I h 7^ for some 1 < h < r — 1. For i = 1, . . . , r let := X^e/iU-uXr 
Then u h ^ and z = ££ =1 u' so z G 

We proceed with the claim about vertices. First we claim that the set of vertices 
of v G W^S™ is contained in V\. Consider any v G W^S™ \ V\. By what we just 
proved, v G" Vo and so 1 < t>j < r — 1 for some i. Since v G" Vi we have that t> + lj 
and v — li are not in Vo hence in the image JV^S", and t> = | ((t> + 1$) + (t> — lj)) 
is a convex combination of them and hence not a vertex. 

Next, consider any v G Vj.. Let i be the unique index where 1 < t>j < r — 1, and 
for t — 0, 1, . . . , r let v l be the vector defined by 



t, if j = i. 



Define a vector h G Z d by 




Then /iw < hv for all m G V, with equality if and only if u G U := {t>°, u 1 , . . . , f r }. 
Therefore W^S*™ D U — {v 1 , . . . ,f r_1 } spans a 1-dimensional face of conv(W / |5'") 
with vertices which are therefore also vertices conv(W T J : S'™). Thus, the ver- 

tices consists precisely of all vectors v G V% having one entry 1 or r — 1. Therefore 
the number of vertices of conv(Wj ! 3'™) is u^(d) = d2 d as claimed. □ 



Combining Lemma 4.2 and Lemma 4.5 we obtain at once the following theorem. 
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Theorem 4.6 The maximum number of vertices of the projection polytope conv(WS) 
of any matroid S C {0, l} n by any {0, l}-valued d x n criteria matrix W satisfies 



d2 d < m(d) < 2^ fa 3 * 3) ) = 0(3^) . 

fc=0 ^ ' 



For d < 2 the lower and upper bounds coincide and we have m(l) = 2 and m(2) = 8. 

Note that already for d = 3 and d = 4 the bounds are far apart, 24 < m(3) < 158 
and 64 < m(4) < 19840, and we do not know the exact values. The value m(2) = 8 
and Theorem 4.1 give the following corollary mentioned in the introduction. 

Corollary 4.7 The maximum number of vertices of conv(WS) for any n, any ma- 
troid S C {0, l} n and any binary 2 x n matrix W , is 8; and max{f(Wx) : x G S} 
can be solved for any convex f by greedily solving 8 linear counterparts over S . 

We conclude this section with an explicit description of the matroid and zono- 
tope realizing the matching lower and upper bounds on m(2) = 8, and the linear 
counterparts used in solving any convex bicriteria matroid optimization problem. 

Example 4.8 Let d = 2, k = 2, r = 3, and n = 8, let Sf be the uniform matroid 
SI = {(1, 1, 1, 0, 0, 0, 0, 0), ... , (0, 0, 0, 0, 0, 1, 1, 1)}, and let Wf be the matrix 



2 



1 1 1 1 
1 1 1 1 



Then, as predicted by Lemma 4.5, conv(W / |S'f) is an octagon, with vertex set 

U := {(1, 0), (2, 0), (3, 1), (3, 2), (2, 3), (1, 3), (0, 2), (0, 1)} , 

homothetic to the zonotope Z(2) = zone{— 1, 0, l} 2 , confirming m(2) = z(2) = 8. 
A set of 8 vectors H which contains, for every vertex u E U, a linear function h u 
which is maximized over conv(Wf 5*|) and Z{2) uniquely at u, is provided by 

H = {(-1, -2), (1, -2), (2, -1), (2, 1), (1, 2), (-1, 2), (-2, 1), (-2, -1)} . 

The polytope conv(W / 2 2 S'f ) and zonotope Z(2), along with the set H, are depicted in 
Figure [2j Given now any n, any matroid S C {0, 1}™ presented by an independence 
oracle, any binary 2 x n matrix W, and any convex function / : IR 2 — > K presented by 
a comparison oracle, to solve ma.x{f(Wx) : x G S} using the algorithm of Theorem 



2.3, proceed as follows: for each of the 8 vectors h G H use the greedy algorithm to 



find Xh G S maximizing hW over S; output that Xh attaining best value f(Wxh)- 
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Figure 2: The octagon conv(Wf 5f) and homothety of zonotope Z(2) 

5 Convex integer programming 

We now turn to convex integer programming, where the set of feasible points is 

S := {x G Z n , Ax = b , I < x < u} . (5) 

We need to define the Graver basis of an integer matrix, see [15] for more details. 
Define a partial order C on IR n where x C y if Xiyi > (that is, x and y lie 
in the same orthant) and \xi\ < \yi\ for i = l,...,n. The Graver basis of an 
integer matrix A is defined to be the finite set Q(A) C Z n of □-minimal elements in 
{x G 7h n : Ax = 0, 0}. For example, the Graver basis of A := (1 2 1) is the set 

0(A) = ±{(2, -1,0), (0,-1,2), (1,0,-1), (1,-1,1)} . 

Let 

110(4)11! := maxilla : x G Q(A)} . 

We have the following Lemma. 
Lemma 5.1 Let S be as in p|) above. Then \\Q(A)\\i is an upper bound on e(S). 
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Proof. Let S be as in (J5]). It has been shown in [5] that the Graver basis of A is 
a set of all edge-directions of conv(.S'); here is the proof, included for completeness. 
Consider any edge e = [x, z] of conv(S') and let h := z — x. Then Ah = and hence 
h = Yli hi is a sum °f some elements hi G all lying in the same orthant, see 

[T5l Chapter 3]. We claim that x + hi G conv(S') for all z. Indeed, /i, G implies 
A(x + hi) = b, and I < x,x + h < u and hi C /i imply I < x + hi < u. 

Now let w G Z n be uniquely maximized over conv(S') at the edge e. Then 
tu/ii = u>(x + /ij) — wx < for all i. But u>/ij = wh = wz — vox = 0, implying that 
in fact whi = and hence x + hi G e for all i. This implies that Ttj is a direction of 
e (in fact, all hi are the same and h is a multiple of some Graver basis element). 

So any edge of conv(S') is parallel to some v G G{A) and hence e(S) < \\G(A) ||i. □ 



In the next two subsections we discuss two important classes of integer programs 
and their applications, where the multicriteria convex problem reduces to a small or 
constant number of linear counterparts, and can be solved altogether in polynomial 



time, using Theorem 2.3 together with Lemma 5.1 and some more ingredients. 



5.1 Convex totally unimodular integer programs 

We have the following theorem on convex totally unimodular integer programming. 



Theorem 5.2 Fix any d. Consider the convex integer programming problem 

max {f(Wx) : x G S} , S := {x G Z n , Ax = b , I < x < u) , (6) 

with A a totally unimodular m x n matrix, b G Z m ; I, u G Z n , W an integer d x n 
matrix, and f : IR d — » K a convex function. Then coiyv(WS) has polynomially many 
O ((n • || W^Hoo)"^ 1 ^) vertices. Further, all vertices of coiav(WS) can be enumerated 
and the convex integer program solved in time polynomial in H^Hoo an d (b,l,u). 



Proof. A circuit of an integer matrix A is a primitive linear dependency on the 
columns of A, that is, a nonzero integer vector satisfying Ax = whose support 
supp(x) is minimal under inclusion and whose entries are relatively prime integers. It 
then follows from Cramer's rule that every circuit x of A satisfies ||x||i < (r+l)A(A), 
where r = rank(A) < min{m, n} is the rank of A and A(A) denotes the maximum 
absolute value of a determinant of a square submatrix of A. 

Now, for a totally unimodular A, the Graver basis G{A) consists precisely of the 
circuits of A, see [151 Lemma 3.19], and A(A) = 1. Therefore for totally unimodular 
A we obtain ||(?(^4)||i < n + 1 and hence by Lemma 5.1 also e(S) < n+1. Moreover, 
if A is totally unimodular then we can realize in polynomial time a linear optimiza- 
tion oracle for S since we can optimize linear functions over S in polynomial time 



17 



using linear programming. The theorem therefore now follows from Theorem |2.3| □ 



A specially appealing application of Theorem 5.2 is to convex multicriteria integer 
transshipment and in particular network flow problems. In the integer transshipment 
problem, the feasible set is S as in ([6]), where A is the totally unimodular vertex-arc 
incidence matrix of the directed graph G = (V, E) underlying the problem; b G is 
the demand vector, with v a supplier if b v is negative and a consumer if b v is positive; 
l e and u e are lower and upper capacities on flow on arc e; W is the criteria matrix 
with WiX the value of transshipment x G Z E under criterion i; and / compromises 



these criteria. We have the following immediate corollary of Theorem 5.2 



Corollary 5.3 Fix any d. Then for every directed graph G, demand b, lower and 
upper capacities l,u, criteria matrix W , and convex function f presented by a com- 
parison oracle, the convex d-criteria integer transshipment problem |#|) can be solved 
in time polynomial in ||W||oo and the binary size (b,l,u) of the rest of the data. 



Theorem 5.2 can be obtained also through the results of [3]. However, therein, 



one needs to repeatedly solve extended linear programs of the form 

{x G lR n , Ax = b , Wx = u, I < x <u} 

over the so called fibers of various u G Z d ; so the totally unimodular or specific 
combinatorial structure of the matrix A is lost. For instance, in solving the convex 
transshipment problem in Corollary |5.3 using our method here, the linear optimiza- 



tion oracle is very efficiently realizable by suitable network flow algorithms, whereas 
the derivation through [1] requires the use of a generic linear programming solver. 

We conclude this subsection with an application to the vector partition problem 
studied in [2J H3]- In this problem we need to partition n given items among p 
players so that player i gets bi items. We are given p x n integer utility matrices 
W l , . . . , W d , where W^j is the utility of item j to player i under criterion k. A 
partition is naturally represented by a p x n matrix X with = 1 if item j is 
assigned to player i and X it j = otherwise. Let W k -X := . Wj^-Xij be the total 
utility under criterion k of partition X. The goal is to find a partition X attaining 
maximum value f{W l ■ X, . . . , W d ■ X) where / : lR d — > R is a convex function. 

Since the partitions form the feasible points of the totally unimodular system 

S = {X G Z p + Xn : J2 X ^ = E X m = ^} , 



we obtain the following immediate corollary of Theorem 5.2 



Corollary 5.4 Fix any d. Then the partition problem with any number n of items, 
number p of players, and b G V with p < n = ^2h> an d px n matrices W 1 , . . . , W d 
and convex function f , can be solved in time polynomial in n and || W^ 1 ^, . . . , || W d ||oo 
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We emphasize that here the number p of players can be variable, whereas in 
previous work (see [13] and references therein) it was assumed to be constant. 

5.2 Convex multi-index transportation problems 

Our result on multi-index transportation problems and other applications follow 
from a more general result about n-fold integer programming. We begin this section 
with a few facts about this theory; see QUI [HJ US] for more details. 

Linear n-fold integer programming is the following problem in dimension nt, 

max {wx : A (n) x = b , I < x < u , x G Z nt } , (7) 

where 





(A, 


Ax ■ 


• A 1 


\ 






A 2 


• 


■ 






A (n) ._ 





A 2 ■ 


• 




(8) 




I 


■ 


■ A 2 


) 





is an (r + ns) x nt matrix which is the n-fold product of a fixed (r, s) x t bimatrix A, 
that is, of a matrix A consisting of two blocks A\, A 2 , with A\ its r x t submatrix 
consisting of the first r rows and A 2 its s x t submatrix consisting of the last s rows. 

Let A be a fixed integer (r, s) x t bimatrix. For any n we write each vector 
x G Z n * as a tuple x = (x , . . . , x n ) of n bricks x 1 G Z*. It has been shown in [I], 
|17j . and [12], in increasing generality, that for every bimatrix A, the number of 
nonzero bricks appearing in any element in the Graver basis Q(A^ n ') for any n is 
bounded by a constant independent of n. So we can make the following definition. 

Definition 5.5 The Graver complexity of an integer bimatrix A is defined to be the 
largest number g(A) of nonzero bricks g l in any element g G Q (A^) for any n. 

We then have the following lemma concerning Graver bases of n-fold products. 

Lemma 5.6 For any integer bimatrix A and all n we have \\G(A^)\\i < g(A)\\G(A 2 )\\ 1 . 

Proof. Consider any Graver basis element g G Q (A^) for any n. Then A^g = 
and hence YH=i A\g % = and A 2 g l = for all i. Therefore (see [151 Chapter 3]), for 
each i, g l can be written as the sum g l = some elements /i M G Q(A 2 ) 

all lying in the same orthant. Let m := k\ + ■ ■ ■ + k n and let h be the vector 

h := (h 1 ' 1 ,...,h 1 ' k \...,h n ' 1 ,...,h n ' kn ) G Z mt . 

Then Y^ij Aih l ' j = and A 2 h hj = for all i,j and hence A^h = 0. We claim that 
moreover, h G Q (A^™ 1 ^. Suppose indirectly this is not the case. Then there is an 
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h G Q (A^) with h C h. But then the vector g G 7L ni defined by g l := Ylf=i h 1 '^ f° r 
all i satisfies g C g contradicting g G Q (AW). This proves the claim. Therefore, by 
Definition |5.5| the number of nonzero bricks /i M of h is at most g(A). So 



Mii = EiiA = E 



1 i=l J=l 



i=l i=l 

Since this holds for all g G G(A^), we get < ^(A)||G(A 2 )||i for all n. □ 

We now have the following theorem on convex n-fold integer programming. 

Theorem 5.7 For every fixed d, p and (r, s) x t bimatrix A, there exists a smallest 
constant v(d,p;A) such that, for every n, every b G 7/ +ns , every l,u G 7L nt , and 
every {0, 1, . . . ,p}-valued d x (n£) matrix W, the projection conv(WS) of the set 

S := {x G Z nt , A^x = 6 , / < x < u) , (9) 

/ias at most v(d,p; A) vertices. Further, all vertices ofconv(WS) can be enumerated, 
and for any convex f : M. d — > R i/ie convex n-fold program meoc{f(Wx) : x E S} 
solved, by solving v(d,p;A) linear counterparts over S, and in polynomial time. 

Proo f. S ince A a nd in particular its second block A2 are fixed, it follows from Lem- 
that S has edge complexity e(S) < \\G(A^) ^ < g(A)\\G(A 2 )\\ 1 



5.1 



and 



5.6 



mas 

bounded by a constant. So it follows from Theorem |2.3| that the number of vertices 
of coiav{WS) is bounded by a constant v(d,p;A) independent of n, b, l,u, W, and 
that the convex problem reduces to solving v(d,p;A) linear counterparts over S. 
Now each linear counterpart is a linear n-fold integer programming problem and is 
solvable in time polynomial in n and {b, I, u), see [H [TTJ, implying the theorem. □ 

The theory of n-fold integer programming has a variety of application, including 
to multi-commodity flows, privacy in databases, and more, see for example [T5J 



Chapters 4-5]. Here we mention only one, but generic, consequence of Theorem 5.7 
to multi-index transportation programs; indeed such integer programs are universal: 
the universality theorem of [6] asserts that any integer program is polynomial time 
liftable to some isomorphic 3 x m x n multi-index transportation program. 

Corollary 5.8 For every fixed d,p,l,m, there exists a constant t(d,p;l,m) such 
that, for every n, every integer line- sums aj^, &i,jfc, Qj, and every {0, 1, . . . ,p} -valued 
d x (7 x m x n) matrix W , the projection comr{WS) of the set of I x m x n tables 

S ■ S X G . ^ %i,j,k ^j,k 1 ^ ^ •Ej,j > k &i,k > ^ , ^i,j,k £ j (10) 

I i 3 k ) 
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has at most t(d,p; I, m) vertices regardless ofn and the line sums. Moreover, all these 
vertices can be enumerated, and for any convex f the problem m&x{f(Wx) : x G S} 
solved, by solving t(d,p; I, m) linear counterparts over S, and in polynomial time. 

Proof. Let A be the (Im, I + m) x Im bimatrix whose first block A\ is the Im x Im 
identity matrix l\ m and second block A2 is the (I + m) x Im vertex-edge incidence 
matrix of the complete bipartite graph K^ m . Then it is not hard to verify that, 
suitably arranging tables in vectors in 'Z lxmxn ^ Z lmn and the line-sums a^, bi^, c% j 
in a right-hand side vector e G 2^ m+n < 7 + m ) = %im+in+mn^ line-sum equations are 
encoded in the system A^ n 'x = e, encoding the multi-index transportation program 
in a suitable n-fold program. The corollary now follows from Theorem 5.7 □ 



The determination of the constants v(d,p;A) in Theorem 5.7 for bimatrices A 



in various applications, and in particular the constants t(d,p;l,m) in Corollary 5.8 



is extremely difficult. Let us briefly discuss the case of binary criteria matrices W 
and define t(d; l,m) := t(d, 1; l,m) to be the largest possible number of vertices of 
conv (WS) for any n, any {0, 1}- valued d x (7 x m x n) matrix W, and any set S of 
/ x in x n tables as in (10). We will bound it from above in terms of the relevant 



Graver complexity. Now, the relevant bimatrix (see proof of Corollary 5.8) is A with 
first block A\ = l[ m and second block A 2 the (I + m) x Im vertex-edge incidence 
matrix of the complete bipartite graph K^ m . Let g(l,m) := g(A) be the Graver 
complexity of this matrix. We note that already g(l,m) is not known for all values 
with 3 < / < m, but bounds are available, [H]. We have the following upper bound. 



Proposition 5.9 The maximum number of vertices of the projection conv(WS) of 
any set oflxmxn tables S in (T(fy by any binary dx (Ixmxn) matrix W satisfies 



d-l 



t(d;l,m) < 2^ 



k=0 



m) min{Z, m} + l) d 
k 



3) 



O ((2g(l,m) min^m})^- 1 )) 



Proof. The second block A 2 of A is totally unimodular so its Graver basis Q(A 2 ) 



consist precisely of its circuits, see proof of Theorem 5.2 and the circuits of A 2 are 
vectors supported on circuits of K^ m with alternating values of 1 and —1. So each 
h G Q{A 2 ) consists of the same number, at most min{/,m}, of —1 and 1 values. 

Now, the proofs of Theorem 2^ and Lemma |5.1| imply that every edge of 
conv(W / 5') is parallel to Wg for some element g G Q(A^). Consider any k = 1, . . . , d 
and let Wk be the corresponding row of W. Write g and Wk as tuples g = (g 1 , . . . , g n ) 
and W k = (Wfc, . . . , Wj}) of bricks with each g i and W{ in U 



Ixm 



g l can be written as a sum g % - 
orthant, with Yli=i ^ — m ) 



Now, for each i, 

Yl'jLi of elements G Q(A 2 ), all in the same 



see proof of Lemma 5.6). Since W is binary we get 



\W k g\ 



n ki 

i=i j=i 



< 



^2^2 \ Wkh l ' j \ < y^min{/, m} < g(l, m) min{Z, m} 



=1 j=i 



=1 j=i 
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So conv(WS) is refined by a zonotope with | ^(2g(7, m) min{Z, m} + 1) — 1 
generators, consisting of one element of each antipodal pair of nonzero vectors in 

{0,±1,±2, . . . ,±g(l,m)mm{l,m}} d . 



Therefore the claimed upper bound now follows from the bound in Lemma 2.2 □ 



Of particular interest is the case of 3 x 3 x n tables, that is, I = m = 3 and 



arbitrary n, which is the smallest where the line-sum equations defining S in (10) 
are already not totally unimodular. Here the Graver complexity is known to be 
g(3, 3) = 9. For bicriteria optimization, that is, d = 2, the bound on the largest 
possible number of vertices of conv(WS) for any n, any {0, l}-valued 2 x (3 x 3 x n) 



matrix W, and any set S of 3 x 3 x n tables with any line-sums as in (10), is 



t(2; 3, 3) < 3024. It would be interesting to determine the value of t(2; 3, 3) exactly. 



6 Open problems 

This article raises many open problems, in particular of combinatorial-geometric 
flavor concerning the determination of the various bounds on the number of vertices 
of projections coiav{WS) of various sets S under small valued and binary W. 

First, as discussed in Section 4, we know very little on the values m(d,p). While 
the relation u(d) < m(d) < z(d) on the maximum numbers m(d) = m(d, 1) and 
u(d) of vertices of the projection of any matroid and any uniform matroid S, re- 
spectively, under any binary W, and the number z(d) of vertices of the zonotope 
Z(d) = zone{ — 1,0, l} d , holds with equality for d = 1 trivially and for d = 2 with 



u(2) = m(2) = z(2) = 8 by Theorem 4.6, we know nothing about higher d. The 



determination of z(d) should certainly be easier and in particular could be obtained 
by brute force computation for small values of d. It would be particularly interesting 
to determine u(d) for all d, while m(d) may be still way out of reach. 

Turning to Section 5.1, an interesting open question related to integer programs 
with totally unimodular defining systems concerns projections of the so-called as- 
signment polytope (often also termed Birkhoff polytope or bistochastic polytope) which 
is the convex hull conv(S') of the set S of n x n permutation matrices given by 

S = {Xe 7L n * n ■ E x « = 1 ' Yl x ^ = ^ ■ 

i j 

Let ad{ri) be the maximum number of vertices of the projection conv(W^S') of the 
assignment polytope by any {0, l}-valued d x (n x n) matrix W into M d . The 
assignment polytope has an exponential number | ^2^=2 (fc) 2 ^'(^ ~ -0' — n^) °^ 



edge-directions, but by Theorem 5.2, for any fixed d, we have that ad{n) is bounded 



by a polynomial in n. It is interesting to determine a^in) and in particular ci2(n). 
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Finally, as discussed in Section 5.2, we know almost nothing even about the 
maximum number t(2;l,m) of vertices of the projection conv(WS) into the plane 
of any set of I x m x n tables of any length n and any line-sums under the pro- 
jection by any binary 2 x (I x m x n) matrix W, except for the upper bound 
(2g(l, m) min{Z, m} + l) 2 — 1 from Proposition 5.9 in terms of the Graver complexity 
g(l, m). The determination of g(l, m) for all 3 < / < m is also open and challenging, 
with the smallest unknown value satisfying g(3,4) > 27. Last, the value t(2; 3, 3) 
for binary projections to the plane of 3 x 3 x n tables is particularly intriguing. 
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